Phonon effects in molecular transistors 
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A rate equation formalism is used to determine the effect of electron-phonon coupling on the 
conductance of a molecule. Interplay between the phonon-induced renormalization of the density 
of states on the quantum dot and the phonon-induced renormalization of the dot-lead coupling is 
found to be important. Whether or not the phonons are able to equilibrate in a time rapid compared 
to the transit time of an electron through the dot is found to affect the conductance. Observable 
signatures of phonon equilibration are presented. 



In recent years it has become possible to fabricate 
devices in which the active element is a very small or- 
ganic moleculei. Such a device may be thought of as a 
'quantum dot': a structure weakly coupled to the macro- 
scopic charge reservoirs ('leads') and small enough that 
the quantization of the energly levels on the dot is im- 
portant. Quantum dots fabricated using conventional 
semiconductor technology have been extensively studied 
experimentallji^ and theoreticallj*^. However, the use of 
small molecules may lead to new physics. In particular, 
as electrons are added or removed from a small molecule 
the shape of the molecule is altered. The energies associ- 
ated with this shape change are not small, and the time 
scales may be comparable to those related to the flow of 
electrons into and out of the molecule. Interesting recent 
data indicate that these effects may lead to interesting 
structures in the conductance spectra of the do f^i^i^ . 

The shape change may be thought of as a cou- 
pling of electrons on the molecule to phonon modes 
of the molecule. The subject of electron-phonon cou- 
pling in quantum dots has received some theoretical 
attentioniSiSiiSiiiiiS, however the existing literature ap- 
pears contradictory, and many papers omit important 
physics or make apparently incorrect statements. For 
example, the phonon-induced renormalization of the 
coupling between the molecule and the leads is of- 
ten neglectediS'ii'i^. Further, some papers assert (in- 
correctly, we believe) that phonon sidebands may be 
observed even in the linear-response conductanca^iiS. 
There is apparently no comprehensive discussion of the 
relation between model parameters and the phonon- 
induced features of the conductance. Further, a cou- 
pling of the dot to leads implies a mechanical coupling 
between the dot and the external world which may al- 
low the phonons on the dot to relax to an equilibrium 
determined by the leads, instead of following the elec- 
trons. The effect of this physics on the conductance has 
apparently not been examined. 

In this paper we present a theoretical treatment of 
a simple model of electron-phonon effects in a molec- 
ular quantum dot which addresses these issues. Our 
treatment is valid in the limit (likely of greatest imme- 
diate experimental relevance) in which the temperature 
is higher than the level broadening induced by coupling 
the molecule to the leads, and is based on a rate equa- 
tion approach with coefficients derived via a golden-rule 



analysis. We present results for average current as func- 
tions of gate and source-drain voltage for two limiting 
cases: phonons uncoupled from the external world and 
responding only to on-dot electrons (phonon equilibra- 
tion slow compared to dwell time of electron on dot), 
and phonons equilibrated to the external world indepen- 
dent of the electron occupation (phonon equilibration fast 
compared to dwell time of electron on molecule). Sub- 
sequent papersi^ will analyze the low-T limit and the 
equilibrium/nonequilibrium crossover. 

For simplicity we consider here the case of a molecule 
with a single relevant level, coupled to a single (Einstein) 
phonon mode and to two leads, which we label as 'left' 
and 'right', described by the Hamiltonian 



H = H, 



dot 



phonon 



Hrr 



(1) 



with 



Hdot = end + Und [ud - 1) + Xujph {h^ + h) nd(2) 

Hphonon — i^phb^b (3) 

Hrm.,el - ^ K ^ + i?.C.) (4) 

a=L,R per 

Here the number of electrons on the molecule, Ud, is 
given by 



nd 



(5) 
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FIG. I: Energy level diagram. fiL,R represent the chemical 
potential in the left and right leads respectively, while e' rep- 
resents the ground state of the singly occupied dot and the 
dashed lines indicate phonon excitations of the singly occu- 
pied dot. The solid line indicates the energy e" + (7 of the 
doubly occupied dot. In this paper we take the U —> oo limit. 
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the parameter U is the charging energy of the molecule, 
the parameter e gives the energy of putting one electron 
into the dot, the V represent the hybridization of the 
molecule and the leads, while the left and right leads are 
characterized by chemical potentials and iir and den- 
sities of states Nl.r- We have defined the zero-phonon 
state to be the ground state when Ud = 0, and we ne- 
glect anharmonicity in the lattice part of the Hamiltonian 
(such anharmonicity is of course induced by the electron- 
phonon coupling and an intrinsic anharmonicity could 
easily be added). 

Fig. 1 shows the energy level diagram corresponding to 
the problem we consider. We choose as zero of energy the 
value e' corresponding to the ground state of one electron 
on an isolated dot, and we consider the behaviour as the 
chemical potentials of the left and right leads are var- 
ied. We specialize here to the large U limit most likely 
to be relevant for small molecules, so restrict the dot 
occupancy to zero or 1. In conventional language the 
source-drain voltage Vsd = , while the gate volt- 

age Vg = ^^^^^ ■ One should include in the model a cou- 
pling between the phonons on the molecule and phonons 
in the leads or elsewhere in the device. Such coupling 
would provide a mechanism for the molecular phonons 
to equilibrate to the leads, and should be included in a 
fully general treatment. Here we study only the limit 
of uncoupled phonons (described by H above) and the 
limit of fully equilibrated phonons, which we obtain as 
described below. 

We shall be interested in the current / flowing between 
the leads, 



dt i—i-pa 



(6) 



as a function of the two chemical potentials. 

The term Hmix leads to broadening of the electron lev- 
els. The resulting widths due to coupling to the left and 
right leads are Vl,r — Nl.rVIj^. We specialize here to 
the case in which temperature is large compared to the 
level broadening: T^^n << T. In this case, the behavior 
may be described by rate equations describing transitions 
between different states of the systemic. A state of the 
system is specified by a number of electrons on the dot, 
and a number of excited phonons, and transitions involve 
changing the number of electrons or phonons or both. 
Before writing the rate equations, it is convenient to per- 
form a standardi^ canonical transformation to eliminate 
the explicit electron-phonon coupling in Eq. ^ Defining 
S = X {J2a '^ct'^ct) (^^ ^ ^) ^i^*! transforming all opera- 
tors via Oe^^ leads to a transformed Hamiltonian 



H' - H'^ot - 
H'dot 



H'„iix with 



= e'ud + LOphO^a + Und {ud - 1) 



(7) 



mix 



L,B. 



aEpa {^Cip^d, + H.C.) (8) 

where the transformed phonon operator a = b—X^^ d\d 
(J, so the phonon ground state depends on the dot occu- 

2 

pancy. e' = e — -^-^ is the 'polaron shift' in the energy for 



adding one electron to the molecule and the interaction 
parameter U is also renormalized, but as we shall focus 
here on [/ — s- cx) we do not write the renormalization ex- 
plicitly here. The crucial phonon renormalization of the 
electron-lead coupling is given by 



X — exp [—A [a) — a)] 



(9) 



After the transformation the state of the system is 
specified by the dot occupancy n and the number of 
phonons q excited above the ground state corresponding 
to the given occupancy n and we shall be concerned with 
the probability that the system has n electrons and 
q excited phonons. Transitions between states are deter- 
mined by applying the golden rule to H'^^^ . The resulting 
transition rate depends on the change in phonon occupa- 
tion number (it also depends on the electron distribution 
in the leads but we write this dependence separately). 

The transition rate involving hopping an electron from 
the dot to lead a and changing the phonon occupancy 
from q (measured relative to the ground state of H'^^^ 
with occupancy n) to q' (measured relative to the ground 
state of with occupancy n — 1) is equal to the tran- 
sition rate involving hopping an electron from the lead a 
to the dot and changing the phonon occupancy from q 
(measured relative to the ground state of H'^^^^ with oc- 
cupancy n — 1) to q' (measured relative to the ground 
state of iJ'jot with occupancy n) and isi^ 



T't,„^Ta\<q'\X\q>\ 



(10) 



The matrix element can be computed by standard 
methods^; its absolute value | < > P = 

is symmetric under interchange of q and q' and is 



y2 



^q<q' 



-A2)'=(gV!)l/2;^|,-,'|g-AV2 



{k)\{{q~k)\{k+\q' -q\)\ 



(11) 



As interesting special cases, let us write the several 
lowest operators: 



Xon 

Xn = (1 - A2) e-^V2 
A21 =V2A(l-f)e-^V2 

= f 1 - 2A2 + 4^1 e-^V2 



(12) 
(13) 
(14) 

(15) 



Observe that for certain values of A some of the matrix 
elements vanish. This unusual behavior is an interference 
phenomenon, which is slightly obscured by the notation. 
A state which has q phonons excited above the ground 
state of the system with n = electrons is a superpo- 
sition (with varying sign) of many multiphonon states, 
when viewed in the basis which diagonalizes the n — 1 
electron problem, and therefore the transition described 
by Xqqi is really a superposition of many different tran- 
sitions, which for some values of A may destructively in- 
terfere. The phonon renormalization of the molecule-lead 
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coupling is apparently omitted, or treated in an average 
manner which neglects the q, q' dependent structure, in 
several recent paper o^^'^^i^^ . 

We are now in a position to write rate equations, fol- 
lowing e.g^. If the phonons are uncoupled from the 
leads we have (fa is the fermi function for lead a). 

K = Ea,,' /. ((9 - 9') ^Ph + Un{n 1)) ri^,P^rli6) 
+ (1 - /a {{q' - q) o^ph + Un{n + 1))) Tl^,P^?^'^ 
- (1 - /a {{q - q') LO.H + Un{n - 1))) V^, ^^P^ 
-fa iiq' - q) LOph + Un{n + 1)) T^, _^P," 

The net current flowing into the dot from lead a is 



I" = - 9) + Un{n 4 

-P"+i (1 - fa iiq - q') Wph + Unin + 



l))r?,,'(17) 



where the sum on n is from to idmax — 1), dmax being 
the maximum occupation of the dot. Note that in steady 
state it follows from the rate equation that ~ —I^ 
which equals the total current across the dot. 

The opposite limit, of phonons equihbrated to the leads 
(assumed here to be at the same temperature) may be 
treated in steady state by forcing the probability distribu- 
tions on the right hand side of Eq.^lto have the phonon- 
equilibrium form P^ = P"e-«'^'"'/^(l - e""'""/'^). In the 
U oo limit this ansatz implies that the probabilty P° 
that the dot is empty is given by. 



E-pa -qujph/T f 
a.q.q' q,ql J a,q,q' 



a,q,q' q,q'^ 



^-/^fa. 



q,q' 



pa 

«,9 



p-qi^ph/T f 
' J a,q,q' 

_(18) 
1 - f 



where /^,, = 1 - /a ((g - q')(^ph), fa,q,q' 

while, Pi = 1 - po. 

In general for both equilibrated and unequilibrated 
cases the steady state rate equations (P„ = 0) consti- 
tute an eigensystem which we solve numerically. From 
these solutions we have computed the current. Rep- 
resentative results are shown in Fig. 2 which plots the 
low-T current as a function of Vgd for two gate voltages: 
Vg = = -IJ-R, upper panel) and Vg = ^ i^iR = 0, 
lower panel), for both equilibrated and unequilibrated 
phonons. 

Steps (broadened by T) in the current associated 
with "phonon side-bands" are observed when the 
source-drain voltage passes through an integer multiple 
of the phonon frequency. However, in the opposite 'linear 
response' limit Vsd ^ (not shown), as Vg is varied we 
find just one main step in the / — y-curve, as Vg passes 
through 0, and only very tiny structures (vanishing as 
g-wph/T^ which is the probability of the dot being empty 
with one phonon excited) when is a non-zero multiple 
of the phonon frequency. This result appears to differ 
from that stated in Ref>i£ who find phonon side bands as 
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FIG. 2: I-V for coupling constant A = 1.0 as a function of Vsd 
for ujph ~ 1 and T — 0.05. Upper panel is for Vg — 0.0, while 
lower panel is for Vg = 14d/2, hr = 0. 



— neglected the fact that the phonon side-bands "float" 
i.e., shift with the fermi level as Vg is changed. 

Fig. 2 reveals on first sight an apparently surprising 
result: for symmetric bias (Vg = 0) and for the coupling 
considered, the current is larger for equilibrated phonons 
than for the unequilibrated case, whereas for the strongly 
asymmetric case i^R = 0), the opposite is true. This is 
surprising because one expects that in the unequilibrated 
case the phonons arrange themselves so as to maximize 
the current. To gain more insight into this phenomenon 
we have calculated the dependence of the ratio of currents 
for unequilibrated and equilibrated phonons on the cou- 
pling A for different degrees of bias asymmetry. We find 
that except for = 0.0 (the most asymmetric case) a 
minimum in the ratio occurs for a A ~ 1. This behaviour 
may be traced back to Ea. ll3ll51 which reveal that higher 
order "diagonal" (n phonon- n phonon) matrix elements 
vanish for a A ~ 1. 

The steps in current may be conveniently parametrized 
by the height (or the area, as the width is simply propor- 
tional to T) of the corresponding peaks Gmax m the dif- 
ferential conductance G = dl/dV. Ratios of peak heights 
(or areas) provide a convenient experimental measure 
of whether the phonons are in equilibrium. At low T, 
the equilibrium phonon distribution corresponds to oc- 
cupancy only of the n = phonon state, so the n-th side 
band involves a transition from the phonon to the n 
phonon state. Therefore the ratios of the peak heights 
or areas are controlled by ratios of |X„oP- In particular 
Eg ns . [T^ 1181 implv that if /i^ — —fiR and T <C Wph, 



QO 

^ max 



equil 



2\X, 



A 



2n 



00 I 



2(n!) 



(19) 



Note that if /i^ 7^ /i^, then /i dependent changes in the 
occupation probabilities lead to additional, and not sim- 
ply characterized n dependence. 

Deviations from this pattern imply non-equilibrium 
Vg is swept at Vsd 0. The difference occurs because Ref phonons. As illustration we display in Fig. 3 Gmax values 
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FIG. 3: Ratio of peak heights for equilibrated and unequi- 
librated phonons and two different coupling strengths and 
fiL ~ —tJ-R. The points for A = 0.5 (open symbols) have 
been multiplied by 10. 



(normalized to the zero frequency peak) for equilibrium 
and non-equilibrium phonons and a weak and strong elec- 
tron phonon coupling. One sees that the non-equilibrium 
case the peak heights display a non-systematic depen- 
dence on electron-phonon coupling and peak index, but 
that in general measurements of the n — 1 and n = 2 
peaks reveal the effect clearly. 



It is also of interest to consider how far out of equi- 
librium the phonon distribution may be driven. Fig. 4 
shows the phonon occupation probabilities for weak and 
strong electron-phonon coupling and Vg = Q. One sees 
immediately that the phonon distribution function is far- 
ther from equilibrium for weak couplings than for strong 
couphngs. We associate this effect to the strong A depen- 
dence of operators Xo„, (Ea. ll2|l which allows the system 
at large A to "jump down" from a highly excited state 
to one of low phonon occupancy. The deviation from 
equilibrium is largest for Vg = for similar reasons. 
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FIG. 4: Phonon probability distributions for two differ- 
ent electron-phonon coupling constants calculated for fiL = 



In summary we have investigated a rate equation 
model for quantum dots with strong electron-phonon 
coupling. We find important effects arising from the 
phonon renormalization of the dot-lead matrix element 
and from whether or not the phonons are able to equi- 
librate. Future papers will address bistabilit}iii*iS and 
also the low-T limit in which the rate equation approach 
is not valid. 
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